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FOREWORD 


of computers. it is assumed that the reader is familiar with 620 programming 
terminology. 


Section 1 provides subroutine entry and exit data, and the formats used to describe the 
subroutines in the sections that follow. 


Section 2 contains programmed arithmetic routines; section 3, elementary function 
routines; and section 4, conversion routines. An index of included routines is provided at 


| 

| 

This manual describes standard subroutines that can be used with the Varian 620 family 
| 

! 

| the beginning of each section. 
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GENERAL DESCRIPTION 


SECTION 1 

GENERAL DESCRIPTION 
| 

| 


1.1 SUBROUTINE ENTRY AND EXIT 


If a subroutine requires only one parameter or argument, programmed entry will be made 
by first loading the desired parameter into the A register and then executing a return 
jump to the subroutine. 


Where more than two input parameters are required, the parameter will be entered into 
the program following the return-jump to the subroutine. The following sequence of 
instructions will be used: 


varian data machines on 
SECTION 1 

| 

| 

| 

| 


Location Instruction Remarks 

| 

EP Return jump Return jump to subroutine. 

Pa 2 Parameter Parameters or parameter locations for 

| subroutine. 

| 

| Poke 3 Parameter Parameters or parameter locations for 
subroutine. 

| P+ 4 Parameter Parameters or parameter locations for 

| subroutine. 

| 

| P+on Parameter Parameters or parameter locations for 

| subroutine. 

P + n + 1Normal return Continuation of program. 

| 

| 

| 

| 

| 
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Y ee FORMAT 


Each routine is organized in the following order: 


IDENTIFICATION 


Symbolic title and description 


PURPOSE 
USE 


Calling sequence or operational procedure 
Arguments or parameters 

Space required (decimal) 

Temporary storage requirements (decimal) 
Error returns or error codes 

Input and output formats 

Sense switch settings 

Accuracy 

Cautions to users 

Equipment configuration 

References 


METHOD OF ALGORITHM 


Items which are not applicable to a particular subroutine have been omitted. 
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| Symbolic Title Description Page 
| ABS Absolute value, floating point 2-17 
| (type real) 
| ABS (FORTRAN version) Absolute value, floating point 2-19 
| (type real) 

IABS Absolute value, fixed point 2-18 

(type integer) 

| IABS (FORTRAN version) Absolute value, fixed point 2-20 
| (type integer) 
ISIG Transfer of sign, fixed point 2-21 
| (type integer) 
SIGN Copy sign 2-22 
| XBTD Fixed-point, single-precision | 
| integer, binary-to-decimal | 
| conversion 2-3 | 
| XDAD Fixed-point, double-precision, | 
| add 2-9 | 
| XDCO Fixed-point, double-precision, | 
| 2's complement 2-8 
XDDI Fixed-point, double-precision, 
| divide 2-15 | 
XDIV Fixed-point, single-precision, | 
| divide 2-7 | 
| | 
| | 
| 
| 98 A 9902 044 2.1 | 
| 
———————————— IE a ae ee NE ene ee TOE | 


SECTION 2 


PROGRAMMED ARITHMETIC 


2.1 GENERAL 


This section contains programmed arithmetic routines, separated into distinct packages. 
Each routine is described according to the format presented in section 1. Items which are 


not applicable to the routine have been omitted. 


2.2 INDEX 


The routines included in this section are listed below, alphabeticaHy by symbolic title, 


along with the page number where they appear. 
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Symbolic Title Description Page 


XDMU Fixed-point, double-precision, 

multiply 2-13 
XDSU Fixed-point, double-precision, 

subtract 2-11 
XDTB Fixed-point, single-precision 

integer, decimal-to-binary 

conversion 2-4 
XMUL Fixed-point, single-precision, 

multiply 2-5 
$FAS Floating add or subtract 2-28 
$FMS Separate mantissa 

(same as $FSM) 2-23 
$FSM Separate mantissa (same as $FMS) 2:23 
$HS Floating-point number to integer 

number 2-24 
SNML Normalize 2-25 
$QK Floating add 2-26 
$QL Floating subtract 2-27 
$QM Floating-point multiply 

(same as $QN) 2-29 
SQN Floating-point divide 

(same as $QM) 2-29 
$QS Integer number to floating-point 

number 2-31 
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IDENTIFICATION 

XBTD Fixed-point, single-precision integer, 
binary-to-decimal conversion. 

PURPOSE 


XBTD converts the absolute value of the integer in the A register, module 10,000, to a 
binary-coded decimal integer in the B register. The input is retained in the A register and 
the X register is unchanged. The output range is 0 through 9999 inclusive. 


ees Td 


Successive division of binary integer by 10,, with concatenation of remainders. 


USE 
Calling sequence CALL XBTD 
Arguments or The binary argument is in the A register 
parameters before and after execution. 
i 
Snace required IR words 
i wy Sees ee wNe oS Vw usw | 
| | 
Temporary storage Four words | 
| required 
| 
| Accuracy Exact 
| 7~~ | ae’ Ly aA - s £ zl By -_ +. a . - «= = 
| Cautions to user An input of -2  wiii set overfiow and 
provide a meaningless result. 

| 
| METHOD 
| 
| 
| 
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IDENTIFICATION 

XDTB Fixed-point, single-precision integer, 
decimal-to-binary conversion. 

PURPOSE 


XDTB converts the binary-coded decimal integer in the A register to a binary integer in 
the B register. The input is retained in the A register with the X register unchanged. The 
output range is +0 through + 9999 inclusive. 


USE 
Calling sequence CALL XDTB 
Arguments or The decimal argument is in the A register 
parameters before and after execution. 
Space required 25 words 
Temporary storage Four words 
required 
Accuracy Exact 
Cautions to users Input is not checked for legal bcd codes, 
but is evaluated as: 
D,* 10° +D, * 10° 4D, * 10! 
+D, * 10° 
where D is a four-bit binary number. 
METHOD 


Successive multiplication of digits by powers of 10 with accumulation. 


B = (10D, +D,) 10 +D,) 10+D, 
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IDENTIFICATION 
XMUL Fixed-point single-precision multiply. 


PURPOSE 


XMUL provides the software version of the (optional) hardware multiply instruction. 


varian data machines (eA) 
SECTION 2 | 
| 


USE 

Calling sequence LDB Multiplier 
| LDA Constant 
| CALL XMUL 
| PZE Address of multiplicand 

Normal return. 

Arguments or On entry, a A = constant to be added 
| parameters to product at 2, 
| B = multiplier. 
| On exit, A, B = double-precision 
| product, 
| X is unchanged. 
| Space required 44,, words 

Temporary storage 

required Two words 
Error returns OV is set (1) if the product is greater 
or codes than 2** (NBIT-1)-1. 
Accuracy Exact | 
Cautions to user OV is reset if there is not an error. 
METHOD 


98 A 9902 044 2-5 
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Recursive addition of muitiplicand with shifting. 
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| 
| 
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Fixed-point single-precision divide. 


PURPOSE 


XDIV provides the software version of one (optional) hardware divide instruction. The true 
remainder and quotient are delivered to the A register and B register, respectively. XDIV 
gives the true result for negative numbers. 


USE 


Calling sequence 


Arguments or 
parameters 


Temporary storage 
required 


Error returns 
or codes 


Accuracy 


Cautions to users 


METHOD 


LDA (Qigh dividend) 
LDB (low dividend) 
CALL XDIV 

PZE (address of divisor) 
Normal return. 


On entry, A, B = double-precision 
dividend. 
On exit, A = remainder, 
B = quotient, 
X is unchanged. 
72 words 
Five words 


OV is set (1) if the dividend is not less 
than the divisor. 


Exact 
This routine produces the true quotient and 


remainder, i.e., 2/1 = quotient of 2 and 
remainder of zero. 


Unsigned, non-restoring divide algorithm. 
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IDENTIFICATION 


XDCO Fixed-point double-precision 2’s complement. 


PURPOSE 


XDCO takes the 2’s complement of the double-precision number in the A and B register. 


The X register is unchanged. 


USE 


Calling sequence 

Arguments or 
parameters 

Space required 


Input and output 
formats 


Accuracy 
Cautions to users 


METHOD 


CALL XDCO 


The A register and the B register contain 
the double-precision argument before, and 
the 2’s complement after execution. 


13 words 


Double-precision numbers are 

stored as two successive data words. 
The first contains the sign and 
high-order 15 bits; the second contains 
the low-order 15 bits and is always 
unsigned. 


Exact 


XDCO may set the overflow register. 


The argument is complemented and the low-order bits are tested for a carry condition. 
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IDENTIFICATION 


PURPOSE 


XDAD adds a double-precision number whose high-order address is in the calling 


| 
XDAD Fixed-point double-precision add. 
sequence to the double-precision numbers in the A and B registers. The X register is 


| unchanged. 
| USE 
| 
| Calling sequence CALL XDAD 
| PZE is the address of the double- 
| precision augend. Normal return. 
| Arguments or The A and B register contain 
| parameters the double-precision added before, 
| and the double-precision sum 
| after execution. 
| Space required 21 words 

Temporary storage Two words 

required 

| Error returns or The overflow is set if a double- 
| error codes precision overflow occurs. 
| Input and output Double-precision numbers are stored 
| formats as two successive data words. The first 
| 
| 
| 98 A 9902 044 2.9 
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contains the sign and high-order 15 
bits; the second contains the low-order 
15 bits and is always unsigned. 


Accuracy Exact 


Cautions to users The sign of the low-order words of each 
double-precision argument must be zero 
to generate the proper carry. Overflow 
flip-flop is set on an overflow. 


METHOD 


Low-order words are added first and any carry generated is added to the high-order sum. 
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XDSU Fixed-point double-precision subtract. 


PURPOSE 


XDSU subtracts a double-precision number (subtrahend) whose high-order address is in 
the calling sequence from the double-precision number (minuend) in the A and B 
registers. The X register is unchanged. 


| 
I 
IDENTIFICATION 
| 


USE 
Calling sequence CALL XDSU 
PZE is the address of high-order bits 
of the double-precision minuend. 
Normal return. 
| Arguments or The A and B registers contain the 
parameiers adoubie-precision subtranend before, and 
| the double-precision difference after 
| execution. 
Space required 23 words 
Temporary storage Two words 
| required 
Error returns or The overflow is set if a double-precision 
| error codes overflow occurs. 
} 
| Input and output Double-precision numbers are stored as 
| formats two successive data words. The first 
\ contains the sign and high-order 15 bits; 
the second contains the low-order 15 bits 
and is always unsigned. 
| Accuracy Exact 
| 98 A 9902 044 2-11 
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Cautions to users The sign of the low-order words of each 
double-precision argument must be zero 
to generate the proper carry. Overflow 
flip-flop is set on an overflow. 
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IDENTIFICATION 
XDMU Fixed-point double-precision multiply. 


XDMU multiplies the double-precision number whose high-order address is in the calling 
sequence times the double-precision number in the A and B register. The X register is 


varian data machines “i 
| 


| 
PURPOSE 
| 
| 
| 
| 
t 
| 


unchanged. 
USE 
Calling sequence CALL XDMU 
PZE is the address of the high-order bits 
of the multiplier. 
Normal return. 
Arguments or The A and B registers contain the 
parameters double-precision multiplicand before 
and the double-precision product after 
| execution. 
| Space required 55 words (without hardware multiply/divide 
option). 49 words (with hardware multiply/ 
divide option). 
| Temporary storage Four words 
required 
Input and output Double-precision numbers are stored as 
formats two successive data words. The first | 
| contains the sign and high-order 15 bits; | 
the second contains the low-order 15 bits 
and is always unsigned. | 
Accuracy 2°" taken as a fraction. 
Cautions to users Operands should be normalized to retain 
| precision. Overflow is reset by XDMU. | 
| 
| 
98 A 9902 044 Die 
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Equipment The hardware multiply/divide option may 
configuration be used; or instead of using the hardware 
option, the XDMU routine can be assembled 
to use the software multiply routine 
XMUL. 
METHOD 


Double-precision addition of partial products. 


(A + a)* (B + b) SAB*2° + Ab*2”° + aBe2' 
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IDENTIFICATION 


XDDI Fixed-point double-precision divide. 


XDDI divides the double-precision number in the A and B registers by the double- 


PURPOSE 
precision number whose high-order address is in the calling sequence. The X register is 


error codes 


unchanged. 
USE 
Calling sequence CALL XDD! 
PZE is the address of high-order bits 
of division. 
Normal return. 
Arguments or The A and B registers contain the | 
parameters doubie-precision dividend before, 
and the double-precision quotient after | 
execution. | 
Space required 83 words (without multiply/divide option). | 
77 words (with multiply/divide option). | 
Temporary storage Six words | 
required | 
| 
Error returns or Overflow is true if a divide fault occurs. | 
| 
| 


Input and output Double-precision numbers are stored as 
formats two successive data words. The first 
contains the sign and high-order 15 bits; 
the second contains the low-order 15 bits 
and is always unsigned. 


| 
| 
| 
Accuracy Accuracy !s + gr taken as a 
fraction. : 

| 
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Cautions to -users Overflow is reset by XDDI. The dividend 
must be less than the divisor. 


Equipment The hardware multiply/divide option may 
configuration be used; or instead of using the hardware 
option, the XDDI routine’ can be assembled 

to use the software divide routine XDIV. 


References XDDI uses XDSU and XDCO. 
METHOD 
A +a A+t+a Aeb 
B+b B B 
2-16 98 A 9902 044 
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IDENTIFICATION 


A 


ABS Absolute value, floating point (type reai). 


PURPOSE 


This routine takes the absolute value of the floating-point (real) quantity in the A and B 
registers, returning the result to the A and B registers. The absolute value of a is defined 
as -a if ais negative, and as a if a is not negative. 


USE 
Cailing sequence CALL ABS 
Arguments or Argument is in the A and B registers. 
parameters 
Space required Six words 
Accuracy No loss of information. 
METHOD 


The method is explained by the coding itself: 


ABS ENTRY Return immediately if not negative. 
JAPes ABS 
CPA One’s complement high order word if 
JMP* ABS negative and return. 


| 
Label Op Code Variabie Comments 
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IDENTIFICATION 
IABS Absolute value, fixed-point (type integer). 


PURPOSE 


This routine takes the absolute value of the signed integer in the A register and returns 
the result to the A register. The absolute value of a is defined as -a if the a is negative and 
a if a Is not negative. 


USE 
Calling sequence CALL IABS 


Arguments or The quantity in the A register is the 
parameters argument. There are no other parameters. 


Space required Seven words 
Accuracy No loss of information. 


METHOD 


The method is explained by the subroutine code itself: 
Label Op Code Variable Comments 
IABS ENTRY 


JAP* JABS Return if argument is positive or 
zero. 


CPA lf argument is negative, one’s 
IAR complement and correct to two's 
complement. 


JMP* IABS Return. 
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IDENTIFICATION (FORTRAN Version) 


ABS Absolute value, floating point (type real). 


This routine takes the absolute value of the floating-point (real) quantity whose address 
follows the CALL instruction, returning the result to the A and B registers. The absolute 


| 
PURPOSE 
value of a is defined as -a if a is negative, and as a if a is not negative. 
| 
| 


ee 


USE 
Calling sequence CALL ABS, ARG 
Arguments or ARG is the address of the argument. 
parameters The result is returned in the A and B 
registers. 
Space required 15 words 
Accuracy No loss of information. 


| 

| 
METHOD 
The coding is the same as the non-FORTRAN ABS except for loading the argument. 
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IDENTIFICATION (FORTRAN Version) 
IABS Absolute value, fixed-point (type integer). 


PURPOSE 


This routine takes the absolute value of the signed integer whose address follows the 
CALL instruction, returning the result to the A register. The absolute value of a is defined 
as -a if the a is negative, and a if ais not negative. 


USE 
Calling sequence CALL IABS, ARG 
Arguments or ARG is the address of the argument. 
parameters The result is returned in the A register. 
Space required, 15 words 
Accuracy No loss of information. 
METHOD 


The coding is the same as the non-FORTRAN IABS except for loading the argument. 
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IDENTIFICATION 
ISIG Transfer of sign, fixed-point (type integer). 


This routine applies the sign of the called (second) parameter to the quantity in the 
accumulator (first parameter). The parameters and result are fixed-point quantities. 


PURPOSE 
USE 
Calling sequence CALL ISIG, REF 
Arguments or The first parameter is located in the 
parameters A register. The second parameter is 
located in core, whose address is in REF. 
Space required 24 words, including two working cells 
(temporary storage). | 
Accuracy No loss of information. | 
METHOD 
Uses $SE. 
} 
98 A 9902 044 2-21 
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IDENTIFICATION 
SIGN Copy sign (floating point). 


PURPOSE 


To set sign of floating point number equal to that of argument. 


USE 
Calling sequence CALL SIGN, REF 
Arguments or Floating point number in A and B registers. 
parameters REF is address of argument. 
Space required 18 words 
Temporary storage Two words 
required 
Input and output Floating point format. 
formats or tables 
Accuracy Exact. 
METHOD 


Sets sign equal to that of argument. Output in A and B registers. Uses $SE. 
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IDENTIFICATION 


$FMS 
Note: 
$FMS and $FSM are two names for the same entry point; use one or 
the other. 
PURPOSE 


| $FSM Separate mantissa (floating point). 

| 

| 

: To separate a positive floating point number into characteristic and mantissa. | 
| 
| 
| 
| 
| 


USE 
| 
| Calling sequence CALL $FMS or $FSM 
Arguments or A and B registers contain fioating 
parameters point number. 
| Space -required 14 words 
\ 
Temporary storage One word 
required 
| Input and output Floating point, input A, B contain fixed 
| formats or tables point mantissa. The characteristic is in 
| bits 15 through 8 of the X register for 
16-bit machines and bits 17 through 10 of 
| the X register for 18-bit machines. | 
| 
Accuracy Exact. 
METHOD 


Output in A, B (mantissa) and X (characteristic) registers. See listing supplied with the 
paper tape of the program. 
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IDENTIFICATION 


$HS Floating point number to fixed-point, single-precision 
integer. 


PURPOSE 


To convert a floating point number to a fixed-point, single-precision integer. 


USE 
Calling sequence CALL $HS, STORE 
Arguments or Number in A and B registers. STORE is 
parameters address of memory where the result is to 
be saved. 
Space required 62 words 
Temporary storage One word 
required 
Error returns or lf number greater than 2 ** 15 or less 
codes than 1, it exits with A and B registers set 
to zero. 
Input and cutput Floating point input. Fixed point 
formats or tables integer output. 
Accuracy 15 bits 
METHOD 


Uses $SE. See listing supplied with the paper tape of the program. 


9.94 98 A 9902 044 


ae aaa ev ee, 


PT varian data machines @ 


SECTION 2 
PROGRAMMED ARITHMETIC 
IDENTIFICATION 
SNML Normalize. 

PURPOSE 

To normalize a double-precision number. 

USE 

Calling Sequence CALL $NML 
| Arguments or Number in A and B registers. 
| parameters 
| Space required 29 words 
| Temporary storage Two words 
| required 
| Input and output Fixed-point format 
formats or tables 
Accuracy 22 bits 

METHOD 

Shifts to sign and tests for sign set. Uses XDCO. Output in A and B registers. Flag for sign 

in X register. 
| | 
| 
| 1 
| | 
| | 
| 
: 
| 98 A 9902 044 2.25 | 
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IDENTIFICATION 
$QK Floating-point add. 


PURPOSE 


To add two floating-point numbers. 


USE 
Calling sequence CALL $QK, REF 
Arguments or A and B registers contain first argument. 
parameters REF is address of second argument. Result 
is in A and B registers. 
Space required Four words 
Input and output Floating-point format 
formats or tables 
Accuracy 22 bits 
METHOD 


Algebraically adds two numbers. 
$QK and $QL use common logic $FAS. $FAS determines if it is an arithmetic addition or 


subtraction and proceeds accordingly. $FAS has a special entry linkage and is used solely 
by $QK and $QL. 
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IDENTIFICATION 
$QL Floating-point subtract. 


PURPOSE 


To compute difference of two floating-point numbers. 


varian data machines (a) 
SECTION 2 
USE 

| Calling sequence CALL $QL, REF 
| 
| Arguments or Minuend in A and B registers. REF is address 
parameters of first word of subtrahend. 
| 
| Space required Four words 
: Input and output See floating-point format. | 
formats or tables 
| Accuracy 22 bits 
| METHOD | 
: Uses $QK. | 
| 
98 A 9902 044 2.27 
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IDENTIFICATION 
$FAS Floating-point add or subtract. 


PURPOSE 
To provide common logic for $QK, $QL. It has a special linkage for use by $QK for $QL. 


USE 
Calling sequence Not for general use. 
Space required 147 words 
Accuracy Exact. 
Cautions for user Not for general use. 
Reference $QK, $QL 

METHOD 


See listing supplied with the paper tape of the program. 
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Floating-point multiply or divide. 


To multiply two floating-point numbers. To divide one number by another. 


USE 


Calling sequence 
Arguments or 
parameters 


Space required 


required 


Error returns or 
codes 


Input and output 
formats or tables 


Accuracy 


98 A 9902 044 


CALL $QM, REF for multiply. CALL $QN, 
REF for divide. 


REF is address of multiplier or divisor. 


133 words 


Seven words 


lf divisor = O, A and B registers set to 
zero and overflow on. 


If result is less than 2** (-200,) or 
greater than 2** (+177.), it returns 
with 0 in A and B registers and overflow 


on. 


Floating-point format. Output in A and B 
registers. 


22 bits multiply 
21 bits divide 


2-29 
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Equipment This routine does not require the hardware 
configuration multiply/divide option; it uses XDMU and 
XDDI which can be assembled to use either 
hardware or software multiply and divide. 


METHOD 
Separate the mantissa and use XDMU for multiply or XDDI for divide. Uses $FMS, $SE. 
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IDENTIFICATION 


PURPOSE 


| 

$QS Fixed-point, singie-precision to floating-point conversion. 
To convert a fixed-point integer to a floating-point. 

| 

| 

| 

| 


USE 
Calling sequence CALL $QS, STORE 
Arguments or Argument™in A register. STORE is address 
parameters of memory where result is to be saved. 
Space required 43 words 
| 
Temporary storage Five words | 
required 
Input and cutput Floating-point format output. Fixed-point 
formats or tables integer input. 
Accuracy Exact. 


Formats the absolute number to floating point and adjusts sign according to input. Uses 
$SE. 


| 
METHOD 
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SECTION 3 
ELEMENTARY FUNCTION ROUTINES 


3.1 GENERAL 


This section contains elementary function routines, separated into distinct packages. 
Each routine is described according to the format presented in section 1. Items which are 
not applicable to the routine have been omitted. 


3.2 INDEX 


The routines included in this section are listed below, alphabetically by symbolic title, 
along with the page number where they appear. 


Symbolic Title Description Page 
ALOG Natural log of floating-point number 3-13 
ATAN Arctangent of floating-point number 3-19 
COS Cosine 3-18 
EXP =xponential 3-16 
POLY Single-precision polynomial 3-11 
SIN Sine 3-15 
SQRT Square root 3-17 
XATN Fixed single-precision arctangent 3-10 
XCOS Fixed single-precision cosine 3-9 
XEXN Fixed single exponential, negative argument 3-5 
XEXP Fixed single exponential, positive argument 3-4 
XLOG Fixed single-precision logarithm 3-3 
XSIN Fixed single-precision sine 3-6 
XSQT Fixed single-precision square root (short) 3-7 
$HE Exponentiation of two integers 3-22 
$PE Exponentiation 3-23 
$QE Exponentiation 3-24 
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PURPOSE 


XLOG computes the natural logarithm of 1 + X, where the single-precision quantity X is in 
the A register. If O0XX<1, the result is returned to theA register, otherwise 
an error exit is taken without further action. Input and output are scaled by2 . 


varian data machines (A) 
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IDENTIFICATION 
XLOG Fixed-point single-precision logarithm. 


USE* 
Calling sequence JMPM XLOG 
JMP (error procedure) 
Normal return. 
Arguments or The argument X is placed in A before 
parameters calling XLOG. 
Space required ZO words 
Error returns or Error return if X is negative. 
error codes 
Accuracy Error is less than 2. machine scale. 
Cautions to users Routine XLOG calls subroutine POLY. 
METHOD 


XLOG uses a Chebychev polynomial of the fifth degree. 
*To compute the natural log of any fixed-point fraction, the following method is used, 
based on the relations LOG (x/y) = LOG (x) - LOG (y), and LOG (x)N = N- LOG (x). 
effectively multiplying the number by 2N. 
2. Remove the sign bit and call XLOG. 


3. Subtract N- (LOGe™ ) fromthe result. 


| 

| 

| 

| 

| 

| 

| 

1. Normalize the number by left shifting until the sign bit is set (N shifts) 
| 

| 

| 
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IDENTIFICATION 


XEXP Fixed-point single exponential, positive argument. 


PURPOSE 


XEXP computes the exponential of X, located in the A register: 


eX 0£X<1. e* is scaled 2. The result is placed in the A register. 
(Also see PURPOSE in subroutine XEXN.) 


USE 

Calling sequence JMPM XEXP 
JMP (error return) 

Normal return. 

Arguments or The argument X is located in the 

parameters A register prior to the call. 

Space required 20 words 

Error returns or An error return is taken without the 

error codes other action if the argument is 
negative. 

Accuracy Error is less than 2. of machine 
scale. 

Cautions to users Note relative scale between input and 
output, and that they differ from scales 
relative to the routine XEXN. System 
subroutine XEXN is called by XEXP. 

METHOD 


The exponential is computed by means of a Chebychev polynomial of the fifth degree. 
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IDENTIFICATION 


XEXN Fixed-point single exponential, negative argument. 


x 2” The result is placed in the A register. (Also see PURPOSE in subroutine XEXP.) The 


PURPOSE 

XEXN computes the exponential of X, located in the A register: eX -1<xS0 eX is scaled 

exponential is split into two subroutines, XEXP and XEXN, to increase scaling flexibility. 
| 


USE 
Calling sequence JMPM XEXN 
JMP (error procedure) 
Normal return. 
Arguments or The argument X is located in the A register 
parameters prior to the call. | 
Space required 18 words | 
| | 
Error returns or An error return is taken without other 
error codes action if the argument is negative. 
Accuracy Error is less than 2." of machine scale. 
Cautions to users Note that scaling conventions. differ be- 
tween subroutines XEXN and XEXP. 
METHOD 


The exponential is computed by means of a Chebychev polynomial of the fifth degree. 
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IDENTIFICATION 


XSIN Fixed-point single-precision sine. 


PURPOSE 


XSIN takes the sine of the quantity X in the A register for range-7w<X<7. The input is” 
scaled by2 . The output is returnedtotheAregister, scaled2 . 


USE 
Calling sequence CALL XSIN 
Arguments or The argument X is in the A register. 
parameters 
Space required 30 words 
Accuracy Error is less than 2)" machine scale. 
Cautions to users XSIN requires subroutine POLY. No test 
is made for7<|x| <4. 
METHOD 


Uses a change of variable to y to reduce range from (-,77) to (-/2, 7/2). The change of 
variable is sin x = siny. 


y=|X- al = if 0 
ye[X- I+ if X<O 


The Taylor sine series, truncated to five items, is used for sin y. 
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IDENTIFICATION 


XSQT Fixed single-precision square root (short). 


XSQT takes the unrounded square root of the quantity in the A register if it is non- 
negative. The result is returned to the A register. The A register is unchanged if the input 
is negative. XSQN is recommended instead, unless there is a hardware divide option. 


H 
PURPOSE 
| 
| 


USE 
Calling sequence JMPM XSQT 
JMP (error procedure) 
Normal return.. 
Arguments or The argument is located in the A register 
parameters before execution. 
Space required 61 words 
Temporary storage Six words 
required 
Error returns or Error return if argument is negative. 
error codes 
Accuracy Error is less than 1.5 x 2. machine 
scale. 
METHOD 


X + 1 = 1/2X) +Alim x; = VA 
2X; 


in the form 


Rte LS Ke ae OK 
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i 
| 
| 


| 
| 
Uses Newton-Ralphson formula 
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where 
AX, = 1/2, A- s) 
Xj 


IfXpo = 1 -2'*° (the maximum positive numeric value of a number ina 
16-bit binary representation) then 4X;< 0 for all steps 


It| AX, | 7 ee ee oe given step, there is no need to take another 
step, as would be required if testing differences of successive x-estimates. A maximum of 
four divide operations makes XSQT less attractive than XSQN (only one divide and one 
short-word multiply) unless automatic divide-hardware is present. 
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IDENTIFICATION 


XCOS Fixed-point single-precision cosine. 


PURPOSE 


XCOS takes the cosine of the quantity X in the A register from range-m XS nm. The 
input is scaled by 2 and the output is scaled by 2 . The output is 
returned to the A register. 


ei re eS a Ri 


varian data machines on 
| 


to -7/2, +7/2. Then cos x = sin y, where y = 2/2 - |X|. The Taylor sine series, 
truncated to five terms, is used for sin y. 


USE 
Calling sequence CALL XCOS 
| 
Arguments or The argument X is in the A register. 
parameters 
Space required 18 words 
| Accuracy Error is less than 2° machine scale. 
| | 
Cautions to users XCOS requires subroutine POLY, no test 
is made form >|x|<4. 
METHOD 


Uses a change of variable to y in order to reduce the range of the variable from (-7, +7) 
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IDENTIFICATION 


XATN Fixed-point single-precision arctangent. 


PURPOSE 


XATN takes the arctangent of the quantity X in the A register, where -1<X<1. The 
input and the output is scaled times 2°. 


ov 


USE 
Calling sequence JMPM, XATN 
Arguments or The argument X is in the A.register. 
parameters 
Space required 14 words 
Accuracy Error is less than 2|" machine scale. 
Cautions to users XATN requires system subroutine POLY. 
METHOD 


XATN uses a Chebychev polynomial of seven terms. This polynomial is adequate for an 
18-bit configuration. 
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IDENTIFICATION 


POLY Single-precision polynomial. 


POLY is a resident utility routine intended primarily to support the fixed-point single- 
precision mathematical subroutines requiring the evaluation of a polynomia! in one 
variable of any finite degree. 


USE 


PURPOSE 
| 

| Calling sequence CALL POLY (list of coefficients, format 
| as below): 
| 

Type code 

List of nonzero coefficients of 


| 
| 
degree greater than 1 
| 
| 


o 


Jarn 
aviv 


Coefficient of degree 1 
Coefficient of degree 0 
Normal return 


>~>oano 


Arguments or The type code is either O or 1. Zero 
parameters denotes a polynomial in all powers; one 
denotes a polynomiai in either odd or 
even powers. 


| 
i 
| 
| 
: The list of coefficients of degree 
greater than one is written highest 
| power first, and may be of any number. 
| d) and e) coefficients must be present. 
Use zero to represent an absent term. 
| 
| 


Space required 47 words 


Temporary storage Three words 
required 


Accuracy The accuracy attainable is close to 
unrounded full single-word precision. 
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However, accuracy obtained depends upon 
correct techniques of scaling and may 
depend on mathematical characteristics 
of the polynomial being evaluated. 


Cautions to users No action is taken if an additive 
overflow occurs during computation of 
the polynomial. Certain arbitrary 
combinations of coefficients may sharply 
reduce the accuracy attained. Missing 
interior coefficients of degrees higher 
than 1 must be approximated by small 
nonzero numbers, unless their absence 
is implied by type code =1. 


METHOD 


The polynomial is evaluated in Horner form. For example: 
4 3 2 
C,X +C; X% + C, X° 4+C, X+C, 
is evaluated as: 
((C4 X +C;) X +C,) X+C,) X + Co 
the parameter list taking the forms 0, C4,C,;, C,, 0, C,, Co. The polynomial 
7 5 3 
C,X +C; X +C, X +C€C, X 
is evaluated as: 


CCX Be RO EOI HO RO 


the parameter list taking the form: 1,C,,C;,C,, 0, C,, 0. 
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IDENTIFICATION 
ALOG Natura! log of floating-point number. 


To compute natural log of a floating-point number. 


eee nn ee en ee eC LS Pt SS PE rt eT 


PURPOSE 
USE 
Calling sequence CALL ALOG, REF 
Arguments or REF is address of argument. 
parameters 
Space required 132 words 
| Temporary storage Eight words 
required 
: Error returns or Exits to $ER if argument = 0. 
| codes 
Input and output Floating-point format. Output in A and B 
formats or tables registers. 
| Accuracy 21 bits | 
| METHOD | 
| Log A = Log ,A * Log 92 | 
| , 
Del | 
i = 0 | 
| FP -V2 | 
ade 
Fo 4+f2 ! 
| | 
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A = F’ * 2D where 1I<F' < 2 


C,4 — , are coefficients of series expansion. Uses $ER, $QS, 
$QOK, $OM, XDMU, XDAD, $FMS, $NML, XDDI, XDSU, $SE routines. 
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IDENTIFICATION 
SIN Sine. 
PURPOSE 


Compute sine of radians ‘in floating point. 


USE 
Calling sequence Call SIN, REF 
Arguments or REF is address (direct or indirect) 
parameters of first word of a floating-point 
number. 
Space required 151 words 
Temporary storage Six words 
Input and output Floating-point format 
formats or tables 
Accuracy 21 bits 
METHOD 


First five terms of Taylor series expansion output in A and B registers. Uses $NML, $QM, 
XDMU, XDAD, $SE, $FMS. 


eee A A 5 NP se ee A re AE hi: ee NPR suns nse neeiatefe AA AR Ae ea oe sniesnnhnerners nin manent shsesmnesmse 
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IDENTIFICATION 
EXP Exponential. 


PURPOSE 


To compute e**A. A is floating-point number. 


USE 
Calling sequence CALL EXP, REF 
Arguments or REF is address of argument A. 
parameters 
Space required 224 words 
Temporary storage Nine words 
required 
Input and output Floating-point format 
formats or tables 
Accuracy 21 bits 
METHOD 


Chebychev approximation uses XDMU, $QK, $QL, $QM, SON, $SE. 
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IDENTIFICATION 


SORT Square root. 


PURPOSE 


To compute square root of a floating-point number. 


een A eS A A cA Seer hs SS 


Newton iteration three times. Uses $SE, XDDI, $FMS. 


USE 
Calling sequence CALL SQRT, REF 
Arguments or REF is address of the argument. 
parameters 
Space required 86 words 
Temporary storage Six words 
required 
Error returns or Exits with zero in A, B if argument 
codes negative and sets overflow flip-flop. 
Input and output Floating-point format 
formats or tables 
Accuracy 21 bits 
METHOD 
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IDENTIFICATION 
COs Cosine. 


PURPOSE 


To compute cosine of angle in floating-point radians. 


USE 
Calling sequence CALL COS, REF. 
Arguments or REF is address of first word of 
parameters floating-point number. 
Space required 19 words 
Temporary storage Two words 
required 
Accuracy 21 bits 
METHOD 


Computes Sime of ( 7/2-A). Uses SIN, $QL, $SE. Output in A and B registers. 
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ATAN Arctangent of a floating-point number. 


PURPOSE 


H 
IDENTIFICATION 
Computes arctangent of radians in floating point. 


LetN = IX lor N =1 X/Y |. The arctangent of N is evaluated by dividing the total range 
OKN< 10 into three intervals: (10 tan 7/24), (tan #w/24, 1), (1, 
10 ). I£N<10°, arctanN=N.IfN>10 , arctanN =7/2. 


USE 
Calling sequence CALL ATAN, REF. 
| 
Arguments or REF is address of the floating-point 
parameters argument. 
Space required 184 word> 
Temporary storage Fight words 
roquirad 
Input and output Floating-point format 
formats 
Accuracy 21 bits 
METHOD 


TAN N = C, N+C,N 4+C,N. 


Continued fraction approximations are used in the remaining intervals. 


| 

| 

| 

| 

| 

| 

| 

| 

| The polynomial approximation in the interval (10° , tanw/24) is: 
| 
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interval (tan 7/24, 1) 


and 


TAN. N = (sign of N) ( 7/2) - N'| D-- 


interval (1, 108) 


where 
C, = 0.99999 99207 
C, = 0.33329 66338 
C, = 0.19574 08066 
A, = 0.23882 29612 
A, = 2.4452 05396 
A, = 1.3247 47223 
B, = 3.9435 29798 
B,; = 1.7982 49626 
D, = 0.99999 92083 
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D, = 0.33328 70775 
D, = 0.06355 00089 
E, = 0.59859 98078 
EF; = 0.39535 44718 


Uses $QM, $QL, $ON, $QK, $SE routines. 


varian data machines nl 
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IDENTIFICATION 


$HE Exponentiation of two integers. 


PURPOSE 


To compute |**J. 


USE 
Calling sequence CALL $HE, REF. 
Arguments or | in A register. REF is address of J. 
parameters 
Space required 29 words 
Temporary storage Two words 
required 
Input and output Fixed-point integers 
formats or tables 
Accuracy 15 bits 
METHOD 


Floats | and uses $PE. Uses $SE, $QS, $HS, $PE. 
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IDENTIFICATION 
PURPOSE 


| 
SPE Exponentiation. 
To compute A**]. 


USE 
Calling sequence CALL $PE, REF. 
Arguments or Argument in A and B registers. REF is 
| parameters address of index I. 
| Space required 34 words 
| 
: Temporary storage Five words 
required 
input and output Floating-point format 
formats or tables 
Accuracy 20 bits 


| 

| 

| METHOD 
| Uses $QS, $QE, and $SE. Floats | and goes to A**B ($QE). 
| 
| 
| 
| 
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IDENTIFICATION 
$QE Exponentiation. 


PURPOSE 
To compute A**B. 


USE 
Calling sequence CALL $QE, REF. 
Arguments or Argument A in A and B registers. REF is 
parameters address of argument B. 
Space required 35 words 
Temporary storage Three words 
required 
Input and-output Floating-point format 
formats or tables 
Accuracy 20 bits 
METHOD 


Uses ALOG, EXP, $SE. 
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| 
SECTION 4 
CODE CONVERSION ROUTINES 


4.1 GENERAL 


This section contains code conversion routines which allow the user to convert from one 
character code, usually associated with a particular peripheral device, to the character 
code of a different device. The three conversion routines described in this section are: 

a. EBCDIC to Hollerith conversion 

b. Hollerith to EBCDIC conversion 

c. EBCDIC to ASCII conversion 


The EBCDIC to Hollerith conversion subroutine (SAQ1) converts an 8-bit EBCDIC 
character in the A register to its equivalent 12-bit Hollerith code in the A register. 


The Hollerith to EBCDIC conversion subroutine (SBO1) converts a 12-bit 029 Hollerith 
character in the A register to its equivalent 8-bit EBCDIC character in the A register. 


| 

| 

The EBCDIC to ASCII conversion subroutine (SCO1) converts an 8-bit EBCDIC character 
in the A register to its equivalent 8-bit ASCII code in the A register. If other than 8-bit 
ASCII code is desired, this routine may be easily modified (see SCO1 subroutine 
description). 
The user should note the following characteristics of these subroutines: 

| a. Requires a VDM 620 series computer with a 16-bit word. 

b. Source statements must be assembled with DAS 8A assembler. 


c. The multiply/divide and extended addressing option is not required. 


This subroutine package is referenced by the following VDM Software part numbers: 
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Source Material 92H0206-001 

Object Materiat 92U0206-001 

Assembly Listing 92L0206-001 
4.2 INDEX 


The routines included in this section are listed below, alphabetically by symbolic title, 
along with the page number where they appear. 


Symbolic Title Description 


Page 
SA01 Convert EBCDIC to Hollerith 4-3 
SBO1 Convert Hollerith to EBCDIC 4-5 
SC0l Convert EBCDIC to ASCII 4-6 
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SAO1 Convert EBCDIC to Hollerith 


PURPOSE 


To convert an EBCDIC character in bits 0 through 7 of the A register to IBM 029 Hollerith 


| 

| 

IDENTIFICATION 
| code in bits 0 through 11 of the A register. 
| 

| 


i 
| 
USE 
Calling sequence P-] LDA vaiue to be converted 
P JMPM SAO] | 
| 
| P +1 
| Peak 2 Any instruction : 
| 
| Arguments or On entry, EBCDIC character in bits 0 | 
parameters through 7 of A register. 
! 
\ 
On exit, X Register unchanged 
B Register unchanged 
A Register converted value in bits 
| O through 11, as follows: | 
| 
| 
| 
| 
| Only one exit exists for this subroutine. 
Return is to P + 2 of the calling program. | 
| 
| Space required 166 words 
Temporary storage Two words 
required 

| 
| 
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Cautions to users This subroutine is not reentrant. Every 
EBCDIC character is convertable. That 
is, there is no error condition asso- 
ciated with this subroutine. 
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! 
| 
IDENTIFICATION 


SBol Convert Hollerith to EBCDIC 


To convert an 029 Hollerith character in bits 0 through 11 of the A register to its 


| 

PURPOSE 
corresponding EBCDIC code in bits 0 through 7 of the A register. 
| 

| 


USE. 
Calling sequence P-] LDA value to be converted 
| P JMPM 
| P+ 1 SBO1 
: P+ 2 Any instruction | 
Arguments or On entry, 029 Hollerith character in bits | 
parameters O through 11 of A register. | 
| On exit, X Register unchanged 
B Register unchanged 
A Register converted value in 
| bits O through 7. | 
| Only one exit exists for this subroutine. 
| Return is to P + 2 of the calling program. 
| Space required 182 words 
! Temporary storage Four words 
| required 
| Cautions to user This subroutine is not reentrant. 
| 
| 
| 
| 
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IDENTIFICATION 


SCO1 Convert EBCDIC to ASCII 


PURPOSE 


To convert an 8-bit EBCDIC character in the A register to its equivalent 8-bit ASCII code 


in the A register. 


USE 


Calling sequence 


Arguments or 
parameters 


Space required 


Temporary storage 
required 


Cautions to users 


P-] LDA value 
P JMPM SCO1 
P+ 2 Any Instruction 


On entry, EBCDIC character in bits 0 
through 7 of A register. 


On exit, B register unchanged. ASCII code 
in bits O through 7 of the A register. 


Only one exit for this subroutine. Return 
is to P + 2 of the calling program. 


84 words 

Two words 

This subroutine is not reentrant. Some 
output devices allow only 7-bit ASCII. 
lf other than 8-bit ASCII is desired, 


this subroutine should be modified as 
follows: 
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cist ify tabi T ” incline 
Fither a. Modify table SCT2 to include 


desired codes 


or b. Insert an appropriate mask 
instruction at location 
SC30 + 1. 
| 
| 
| | 
| | 
| 
| | 
| 
| 
| 
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EVALUATION QUESTIONNAIRE 


TITLE 
MANUAL NUMBER 


The purpose of this questionnaire is to provide suggestions about how the manual can be improved when it is revised. 
It is the goal of the Technicai Publications Department to make each manual as useful as possible and at the same 
time eliminate material that is of no practical value to the user or Customer Service Representative in acquiring 
initial knowledge of, and in maintaining, the equipment in the field. You, as the person working most closely with 
the manual and the equipment, can best provide the input needed by the writer to make the best possible manual for 
your use. 


1. Please complete the following chart. 


CHAPTER/ 
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